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ABSTRACT 

We study the three-dimensional distribution of matter at z ~ 2 using high reso- 
lution spectra of QSO pairs and simulated spectra drawn from cosmological hydro- 
dynamical simulations. We present a sample of 15 QSOs, corresponding to 21 baselines 
of angular separations evenly distributed between ~ 1 and 14 arcmin, observed with 
the Ultraviolet and Visual Echelle Spectrograph (UVES) at the European Southern 
Observatory- Very Large Telescope (ESO-VLT). The observed correlation functions of 
the transmitted flux in the Hi Lyman-a forest transverse to and along the line of 
sight are in agreement, implying that the distortions in redshift space due to peculiar 
velocities are relatively small and - within the relatively large error bars - not signif- 
icant. The clustering signal is significant up to velocity separations of ~ 300 km s~^, 
corresponding to about b comoving Mpc. Compatibility at the 2 <j level has been 
found both for the Auto- and Cross-correlation functions and for the set of the Cross 
correlation coefficients. The analysis focuses in particular on two QSO groups of the 
sample, the Sextet and the Triplet. Searching for alignments in the redshift space 
between Lyman-a absorption lines belonging to different lines of sight, it has been 
possible to discover the presence of a wide Hi structures extending over about ten 
Mpc in comoving space, and give constraints on the sizes of two cosmic under-dense 
regions in the intergalactic medium, which have been detected with a 91% and 86% 
significance level, respectively in the Sextet and in the Triplet. 

Key words: intergalactic medium, quasars: absorption lines, cosmology: observa- 
tions, large-scale structure of Universe 



1 INTRODUCTION 

The understanding of the Lyman-a forest has dramatically 
improved in the recent decade, both on the theoretical and 
the observational side. Semi-analytical and hydro- dynamical 
simulations have outlined a new picture where the Lyman- 
a forest is due to the fluctuations of the intermediate and 
low-density intergalactic medium (IGM), arising naturally 
in the hierarchical process of structure formation. Relatively 
simple physical processes impact on the thermal state of the 
gas, which, on scales larger than the Jeans length, effectively 
traces the underlying distribution of dark matter. Support 
for this scenario is given by the satisfactory reproduction 
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by semi-analytical and hydro-simulations of many proper- 
ties of the Lyman-a forest (from the column density and 
the Doppler parameter distribution to the number density 
and effective opacity evolution) derived from the analysis of 
high resolution, high signal-to-noise ratio (SNR) QSO spec- 



tra obtained at 8- 10m class telescopes (e.g., Dave et al. 1999 
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Janknecht et al. 2Q06| . A step forward in the study of the 
IGM with QSO absorption spectra is represented by the use 
of multiple lines of sight at close angular separations, which 
allows information about the transverse direction to be ob- 
tained. Common absorption features observed in the spectra 



of multiply lensed quasars (e.g., Smette et al. 
close quasar pairs (e.g., [D'Odorico et al.||1998 



1992) and of 



Aracil et al. 



2002) have provided evidence that the Lyman-a absorbers 



have dimensions of a few hundred kpc, in agreement with 
the predictions of simulations. Recently, lensed and more 
widely separated QSO pairs have been used to recover the 
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kinematics of the gaseous cosmic web dRauch et aT]|2005| ), 
confirming that the Hubble expansion and gravitational in- 
stability are the main processes influencing the Lyman-a 
forest gas. 

A critical test of the nature of the Lyman-a absorbers, 
as proposed by simulations, comes from the determination 
of their spatial distribution properties. This has been done 
by analysing a great number of uncorrelated QSO lines of 
sight and computing the flux correlation function by averag- 
ing over many spectra (e.g., Tripp et al.|1998 Savaglio et al. 
1999 Croft et al.| 2002t. In this observational approach, how- 
ever, the three-dimensional information is convolved with 
distortions in redshift space, due to peculiar motions and 
thermal broadening. Multiple lines of sight at small angu- 
lar separations offer an invaluable alternative to address the 
spatial distribution of the absorbers, enabling a more direct 
interpretation of the observed correlations. 

The final goal of this work is to investigate the distribu- 
tion properties of matter in the IGM applying the modern 
interpretation of the Lyman-a forest to a sample of close 
QSO groups. The computation and comparison of the flux 
correlation function along and across the lines of sight was 
carried out in a previous paper dP'Odorico eraL][2006l Pa- 
per I). In this paper, we update the previous results using 
higher SNR spectra and, in particular, we investigate the 
coincidences of absorbers among three or more close lines 
of sight in order to detect cosmological structures extending 
to large scales. Finally, we look for extended under-dense 
regions in the IGM using multiple lines of sight. 

Note that tomographic studies based on Lyman-alpha 
lines similar to those presented here will be of great impor- 
tance in the near future due to the large numbers of quasar 
spectra that will be collected using low and medium resolu- 
tion spectrographs (BOSS survey and the X-Shooter instru- 
ment for example) that aim at constraining spatial matter 
correlations in the transverse and longitudinal directions to 
ultimately measure baryonic acoustic oscillations and per- 
form the Alcock-Paczyinski test at high redshift. 

The paper structure is the following: in Section 2 we 
describe the observed data sample, the reduction procedure 
and the simulated spectra. Section 3 is devoted to the com- 
putation of the auto- and cross-correlation functions and 
the cross correlation coefficients of both the observed and 
simulated spectra and to their comparison. In Section 4, we 
analyze the coincidences in redshift space of the H i and C IV 
absorbers in order to detect structures extending over sev- 
eral comoving Mpc, while the under-dense regions common 
to multiple lines of sight are studied in Section 5. The con- 
clusions are then drawn in Section 6. Throughout this paper 
we adopt Qom — 0.26, Oqa — 0.74, and h — Hq/[72 km s~^ 
Mpc-^). 



2 DATA SAMPLE 

The QSO pairs and groups forming our sample are the same 
as Paper I and are described in Table [l] The patterns on the 
sky of the two QSO groups called the Sextet and the Triplet 
are reported in Fig. [l] 

We were allocated 17 hours with UVES at the VLT 
in order to increase the SNR of the spectra of objects Tl 
and T2 of the Triplet and SI, S2, S5 and S6 of the Sextet. 



Table 1. Characteristics of the observed QSO spectra. 





Object 




z 


Mb 


Lyman-o; 
range 


SNR 
per pixel 


Pair A 


PAl 


2.645 


19.11 


2.094-2.585 


8-12 




PA2 


2.610 


19.84 


2.094-2.550 


3.5-6.5 


Triplet 


Tl 


2.041 


18.20 


1.633-1.991 


3-15 




T2 


2.05 


18.30 


1.592-1.999 


4-15 




T3 


2 


053 


18.10 


1.665-2.002 


2.5-7 


Sextet 


SI 


1 


907 


19.66 


1.665-1.859 


2-7 




S2 


2 


387 


19.53 


1.858-2.331 


3-8 




S3 


2 


102 


19.31 


1.633-2.051 


4-12 




S4 


1 


849 


19.59 


1.575-1.802 


3-9 




S5 


2 


121 


18.85 


1.633-2.069 


4-12 




S6 


2 


068 


20.19 


1.592-2.017 


3-10 


Pair U 


UM680 


2 


144 


18.60 


1.653-2.092 


6.5-17 




UM681 


2 


122 


19.10 


1.634-2.070 


7-17 


Pair Q Q2343+12 


2 


549 


17.00 


1.994-2.490 


13-23 




Q2344+12 


2 


773 


17.50 


2.183-2.711 


12-18 



The new observations were reduced with the UVES pipeline 
following 



standard procedure (Ballester et al. 2000). A 
pre-filtering of the cosmic rays of the blue band frames was 
necessary for the faintest spectra to carry out the optimal 
extraction properly. 

Then, the heliocentric and vacuum wavelength cor- 
rections were performed and the new spectra were com- 
bined with the old ones. The final spectra have resolution 
R ^ 45000, while the SNR per pixel varies on average be- 
tween 3 and 12 in the Lyman-a forest and between 6 and 
20 in the Civ forest (see Tab. [l] for details). 

The estimate of the continuum level, in particular in 
the Lyman-a forest region, is a very delicate step in the 
process of spectra reduction. Procedures realized up to now 
in order to determine the continuum position through auto- 
matic algorithms do not give satisfactory results. Following 
the same procedure adopted in Paper I, we fitted the regions 
free from clear absorption with a spline polynomial of third 
order. The limitations introduced by the uncertainty in the 
true continuum level should play a minor role in the compu- 
tation of the cross-correlation function with respect to the 
the case of single line of sight analysis. This is because the 
power on scales of the continuum fluctuations is uncorre- 
lated between the different lines of sight. Thus, probing the 
transverse direction could potentially allow to measure the 
matter distribution at scales larger than those probed along 
the line of sight (see Viel et al.|2002 for a discussion in the 
case of power spectra). 

All the lines falling in the Lyman-a forest have been 
fitted with a Voigt profile via minimization. The lines 
with an equivalent width (EW) lower than three times the 
related EW uncertainty have been removed from the list of 
fitted lines. The metal lines have been identified first looking 
for the most common doublets (e.g. Civ, Si IV, Ovi, and 
Mgii). Then, we have searched for other common transitions 
(e.g. Sim, Sin, Cii, Fell) at the redshift of the previously 
determined systems. 
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Figure 1. Relative positions, U magnitudes and redshifts of the QSOs composing the Sextet (left plot) and the Triplet (right plot) in 
our sample. 



Our sample provides 21 QSO pairs with angular separa- 
tions uniformly distributed between ^ 1 and 14 arcmin, cor- 
responding to comoving spatial separations between ^1.4 
and 21.6 Mpc. The median redshift of the Lyman-a for- 
est is z ^ 1.8. This is the largest sample of high-resolution 
spectra of QSO pairs ever collected, unique both for the 
number density - we have six QSOs in a region of ^ 0.04 
deg^ - and the variety of line of sight separations investi- 
gated. 



2.1 Simulated spectra 

In order to both assess the nature of the Lyman-a forest 
inferred from simulations and to constrain the cosmological 
scenario of the same simulations, we compared the results 
obtained for our sample of observed QSO spectra with anal- 
ogous results for a sample of mock Lyman-a forests. 

The details of the adopted simulations can be found 
in Paper I, here we provide only the basic information. 
We used simulations run with the parallel hydro- dynamical 
(Tree SPH) code GADGET-2 ( jSpringel et"aL]|2001| [Springel 
|2005 | . The simulations were performed with periodic bound- 
ary conditions with an equal number of dark matter and gas 
particles and used the conservative 'entropy-formulation' of 
SPH proposed by Springel & Hernquist (2002). Radiative 



cooling and heating processes were followed for a primor- 
dial mix of hydrogen and helium. We assumed a mean UV 
background produced by quasars and galaxies as given by 



Haardt & Madau (1996) with helium heating rates multi- 



plied by a factor 3.3 in order to fit observational constraints 
on the temperature evolution of the IGM. More details can 



be found in Viel et al. (2004) 



The cosmological model corresponds to a 'fiducial' 
ACDM universe with parameters Qom = 0.26, Ho a = 
0.74, Q q^ = 0.0463 and Ho = 72 km s"^ Mpc"^ (the B2 
series of Viel et al.|2004 ) . We have used 2 x 400^ dark mat- 



Spergel et al. 2003 Viel 



ter and gas particles in a 120 comoving Mpc box. The 
gravitational softening was set to 5 kpc in comoving 
units for all particles. We note that the parameters cho- 
sen here, including the thermal history of the IGM, are in 
perfect agreement with observational constraints including 
recent results on the GMB and other results obtained by the 
Lyman-a forest community (e.g, 
et al.||2004||Seljak et al.|2005 ). 

The z = l.S output of the simulated box was pierced by 
three lines of sight in order to obtain 50 triplets of spectra 
carefully reproducing the observed Triplet mutual separa- 
tions and spectral properties. The same was done for 50 
sextets of lines of sight reproducing the observed Sextet and 
50 pairs of lines of sight at the same angular separation as 
Pair U. 50 different realizations of Pair A spectra and of Pair 
Q were obtained from the output box at redshift z = 2.4. 

Finally, we added to the simulated spectra both instru- 
mental broadening and a Gaussian noise in order to repro- 
duce the observed average SNR (per pixel): SNR = 5 for 
the Triplet, the Sextet and Pair A; SNR = 9 for Pair U and 
SNR = 15 for Pair Q. 



3 CORRELATION FUNCTIONS OF THE IGM 

In this section, we discuss the flux correlations in the ab- 
sorption spectra of our sample of QSO pairs. The statistical 
quantities are the same as those already computed in Paper 
I, here we verify the effect of having an increased SNR. 

On the basis of the interpretation of the Lyman-a for- 
est as due to a continuous density field with a one-to-one 
correspondence between density and transmitted flux, we 
computed the correlation properties of the transmitted flux 
in QSO lines of sight and regarded them as indicators of 
the correlation properties of matter in the IGM. We se- 
lected in each normalized spectrum the region between the 
Lyman-/3 emission (or the shortest observed wavelength. 
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when the Lyman- /3 was not included in the spectrum) and 
5000 km s~^ from the Lyman-a emission (to avoid proxim- 
ity effect due to the QSO). Absorption Unes due to ions of 
elements heavier than hydrogen contaminate the Lyman-a 
forest and can give spurious contributions to the clustering 
signal (see Kim et ar||200 4* for a discussion in the case of 
single lines of sight). We flagged and removed the spectral re- 
gions where metal lines and Lyman-a absorptions of damped 
and sub-damped systems occurred inside the Lyman-a for- 
est. 

Given the normalized transmitted flux, /, as a function 
of the velocity v\\ along the line of sight and the angular 
position on the sky, we define Sf = {f — f), where the 
average fiux, /, is computed for every spectrum as the mean 
of the transmitted fiux over all the considered pixels in that 
spectrum. We neglected the redshift evolution of the average 
transmitted fiux in the Lyman-a forest of the individual 
spectra, which translates into the redshift evolution of the 
mean Hi opacity of the Uni verse (|Kim et al.||2002] [Schaye 



eraL]|2003l |Viel et al.||2004l [Kirkman et al.||2b05| |Fauch~ 



Giguere et al. 2008), because we verified that its effect on 
the correlation function is negligible. By means of this new 
field, we could then define and compute three useful tools 
for the investigation of the Lyman-a correlation properties: 
the Auto and Cross correlation function and the set of the 
Cross correlation coefficients. 



3.1 The Auto-correlation function 

The unnormalized Auto correlation function (Auto CF) of 
the fiux along the line of sight is defined as: 



d{Av^^) = +A«||)>, 



(1) 



following previous studies on the same subject (e.g. McDon- 



ald et al. 2000 Rolhnde et al. 2003 Becker et al. 2004). 



The Auto CF for our sample of QSO spectra was obtained 
by averaging over all the pixels of all the QSOs. The results 
were binned in 50 km s~^ velocity bins. The Auto CF for the 
simulated spectra was computed as the arithmetic mean of 
the correlation functions obtained for 50 realizations of the 
observed sample and the error is the corresponding standard 
deviation. It is important to recognize that the computed er- 
ror bars for the simulated are strongly correlated. This is 
due to the fact that every pixel contributes to the correlation 
function in several velocity bins. 

In Fig. [2] we show the comparison of the observed and 
simulated Auto CF in velocity space. With the improved 
SNR with respect to Paper I, the agreement between the two 
correlation functions has increased, weakening any evidence 
for extra clustering in the regions occupied by QSO pairs. 



3.2 The Cross-correlation function 

In this section we exploit the capabilities of our sample of 
QSO pairs by determining the clustering properties of the 
IGM across the lines of sight. The great advantage with 
respect to the correlation function along the line of sight, 
in particular for a sample like ours showing a large variety 
of pair separations, is that we have the guarantee of sam- 
pling true spatial separations between the pixels, the effect 
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Figure 2. Comparison of the observed and simulated Auto- 
correlation functions: solid dots refer to the observations and 
empty dots to the simulations. Simulated data are slightly shifted 
in velocity for clarity. 



of peculiar velocities being negligible or absent. As a first ap- 
proach, we computed the Cross-correlation function (Cross 
CF) extending in a natural way the procedure adopted for 
the Auto CF. 

Every pixel along the line of sight is considered as an 
element of the density field at the QSO angular position in 
the sky and at a distance from the observer (comoving along 
the line of sight) corresponding to the redshift of the pixel: 



^11 (^) 



^0 Jo 



dz' 
E{z') 



(2) 



where E{z) describes the evolution of the Hubble parameter 
as a function of the redshifd 

In the definition of the comoving distance there is the 
implicit hypothesis that peculiar velocities give a negligible 
contribution to the measured redshift in the Lyman-a forest 
dRauch et al.|2005[ Paper I). The Cross CF of the transmit- 
ted fiux between two lines of sight at angular separation 
is defined as 



^{(Ar) = (5/(^,r||,i)5/(^ + A0,r||,2)), 



(3) 



where, Ar = 



+ 2 ~ 2 r\\^i r\\^2 cos A^) is the spatial 
1 along one line of sight and 



separation between pixel 1 at r 
pixel 2 at r\\^2 along the paired line of sight. 

We also computed the Cross CF for the sample of mock 
spectra. The simulated spectra are characterized by the red- 
shift of the output box and a velocity extent. In order to as- 
sign a redshift value to every pixel, we gave the central pixel 
of every spectrum the redshift of the corresponding output 
box, then we numbered the pixels one by one, transform- 
ing the velocity size of the pixel into a redshift size. Once 
the redshifts were determined pixel by pixel, we followed the 
same procedure adopted for the observed spectra for the 50 
simulated samples and computed the average Cross CF and 



^ E{z) = Y^f^om (1 + ^)3+^^0 
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Figure 3. Comparison of Cross-correlation functions: solid dots 
are related to the observations and empty points to the simula- 
tions. Simulations data are slightly shifted in comoving distance 
for clarity. 



0.020 



-0.00^ 





4 6 

Ar [comoving h"^Mpc] 



Figure 4. Comparison between Auto and Cross-correlation func- 
tions, shown with triangles and circles, respectively. Errorbars 
show the uncertainties related to the simulated samples. 



its 1 a standard deviation. The result of our computation is 
shown in Fig. |3] compared with the observed Cross CF of the 
pairs. The two functions show very good agreement (at the 
1 a level) and a significant clustering signal up to comoving 
separations of ^ 4 comoving Mpc. 

We also compared the Cross CF with the previously 
computed Auto CF. As can be seen in Fig. [4] the two data 
series show a consistency at the 1 a level, supporting the 
hypothesis that the peculiar velocity along the lines of sight 
is negligible. 



3.3 The Cross correlation coefficient 

A measure of the transverse clustering properties of the IGM 
which is less affected by peculiar velocities is the set of the 
flux Cross correlation coefficients (CCC), 

xiiAe) = {Sfie,v^^)5f{e + Ae,v^^)}, 



(4) 



where, every pixel along one line of sight is correlated with 
the one face-to-face in redshift space along the paired line 
of sight and the result is averaged over all the pixels in the 
common redshift interval. 

Every pair of QSOs at angular separation gives one 
value of Xx i^^)^ ^ind a sample with several pairs at different 
separations, as in our sample, gives an estimate of the corre- 
lation function. At a given redshift, the angular separation 

corresponds to a velocity separation A^;^ = cF{z)AO, 
where c denotes the speed of light, and F{z) is a dimension- 
less function of redshift that includes all the dependence on 
the global cosmological metric. In the cosmological model 
that we have adopted: 



F{z) 



Eiz)Qdz'/E{z')] 
(i + z) 



(5) 



where z is the mean redshift of the overlapping Lyman-a re- 
gion from a pair. We computed the range of velocity separa- 
tions covered by each of our pairs of spectra then we grouped 
the pairs in velocity bins of variable width and computed 
the average CCC for every group. Given the small number of 
pairs in each group (a maximum of three QSO pairs) the un- 
certainties associated with these determinations can be only 
computed by applying a simple error propagation. However, 
these uncertainties would account only for the pixel statis- 
tics and the noise associated with the transmitted flux, but 
would be not representative of the true error due to the cos- 
mic variance. In the case of the simulated spectra, we had 
50 realizations of each of our QSO pairs, so we could obtain 
in every velocity interval deflned for the observed pairs an 
average CCC with its error, that in this case is the standard 
deviation of the distribution of values. 

The results are shown in Fig. [5] The CCC computed 
from the observations are in excellent agreement, at the 1 a 
confldence level, with the simulations. 

A comparison can be done with the same results ob- 
tained in Paper I: thanks to the increased SNR for some 
of the spectra in the present sample, in particular the Sl- 
S3 pair, the enhanced clustering signal measured in Paper I 
with the cross correlation coefficient at a transverse velocity 
separation Av± ^ 500 km s~^ is no longer significant. 

We conclude this section by comparing the CCC with 
the Cross CF. Both estimators measure the clustering prop- 
erties of the IGM in the transverse direction with respect to 
the line of sight. However, the Cross CF could be affected 
by the peculiar velocities in the ^-direction while their effect 
should be negligible for the CCC. 

Concerning the CCC, the angular separation A^ be- 
tween two QSO lines of sight was transformed into a comov- 
ing spatial separation, Ar, with the formula 



Ar : 



cAO 



dz' 
E{z') 



(6) 



where E[z) was defined above and here z is the mean redshift 
of the Lyman-a interval considered for each pair. Assigning, 



6 M. Cappetta et al 



0.020 



0.005 



200 



400 



600 800 
A [km s~^] 



1000 



1200 



1400 



Figure 5. Comparison of the Cross correlation coefficients for 
our sample of observed spectra (solid dots) and of simulated 
ones (empty circles) as a function of the velocity separation, 
Av± = cF{z)A6, corresponding to the angular separation, A^, 
of the QSO pairs. Simulated data are slightly shifted in veloc- 
ity for clarity. Error bars on the observed values along the x-axis 
represent the velocity range covered by the considered pairs. 
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Figure 6. Cross correlation function of the transmitted flux, , 
for our sample of pairs and groups of QSOs (empty triangle, 
slightly shifted in Ar) compared with the observed Cross correla- 
tion coefficients, Xx (squares), as a function of spatial separation 
(see text). Error bars both on Xx ^^^^ been determined 

from simulations. 



at both the CCC and the Cross CF observed data points, 
the error bars computed from the simulations, there is an 
agreement within 1 a between the two series of points (see 
Fig.[6|, confirming the hypothesis that the peculiar velocity 
along the lines of sight are negligible. The large variations 
from one data point to the other in the CCC should be due 
mainly to the small number of QSO pairs (between one and 
three) contributing to each point. On the other hand, the 
smoothness of is artificially increased by the fact that 
the values in the different bins are not independent. 



All the results reported in this section are in agreement 
with the same quantities derived from other larger sample 
of QSO pairs (see for e.g. Rollinde et al.|2003 Coppolani et 
|^1.„2QQ6J . 



4 COINCIDENCES OF THE ABSORPTION 
LINES 

4.1 The Lyman-a line coincidences 

The Cross CF and the CCC are two powerful tools which 
allow us to measure the statistical properties of the IGM 
distribution in the real 3D comoving space. By definition, 
they account for the correlation between pairs of Lyman-a 
forests in close QSO spectra. What we want to study now is 
the simultaneous correlation among three or more Lyman-a 
forests. A natural way to extend the Cross CF of the trans- 
mitted flux to more than two lines of sight is the addition of 
a third parameter, related to the third QSO at angular sep- 
aration A^' with respect to the first. This operation leads 
to a function of three variables, the three reciprocal comov- 
ing distances between the three considered pixels, which is 
computation time demanding and difficult to interpret. 

In the previous section, we have taken advantage of the 
interpretation of the Lyman-a forest which allows us to map 
directly the transmitted flux along the line of sight into the 
IGM density field. A boost in the signal of the Cross CF 
between two lines of sight is due to the presence in redshift 
space of two aligned, or very close, Lyman-a lines belonging 
to the two considered spectra. On this basis, we can provide 
a measure of the cross correlation between three Lyman-a 
forests by searching for triplets of Lyman-a lines, belonging 
to three different spectra, aligned in redshift space within a 
given velocity window. 

This kind of analysis has been applied to the Triplet and 
to all the combinations of 3 QSOs that could be formed with 
the Sextet. The adopted procedure has been the following: 

1. The lists of Hi Lyman-a lines compiled for the QSOs 
in our sample were considered in the redshift range between 
the Lyman- /3 emission (or the shortest observed wavelength, 
when the Lyman- /3 was not included in the spectrum) and 
5000 km s~^ from the Lyman-a emission (to avoid proximity 
effect due to the QSO). 

2. Each pair of lines with a velocity separation Av ^ 100 
km s~^ has been replaced by a single line with central wave- 
length equal to the average value of the parent lines weighted 
on the EW. This velocity threshold has been chosen on the 
basis of the characteristic width of Lyman-a lines, ^ 25 — 30 
km/s (see e.g. Kim et al.|2002( ). Furthermore, this is also the 
velocity scale corresponding to the Jeans length, which sets 
the characteristic dimension of H i absorbers. 

3. Triplets of lines, each one belonging to a different line of 
sight, have been considered and the velocity difference be- 
tween the largest and smallest redshift has been computed. 
This operation has been done for the Lyman-a lines in the 
three lines of sight of the Triplet and in all the triplets of 
lines of sight (20 possible combinations) provided by the 
Sextet. Then, all the measures of velocity difference lower 



than 1000 km s~ have been divided into velocity bins of 
100 km s~^ and the related histogram with the number of 
occurrences for each bin has been computed. 
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Figure 7. Three point probability excess (PE3) as a function of 
the velocity difference of the three Lyman-o; lines (dots). Values 
referring to a significancy level of 1, 2 and 3 a levels are reported 
too (dashed, dot-dashed and dot lines respectively). 
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Figure 8. Four point probability excess (PE4) as a function of 
the velocity difference of the three Lyman-o; lines (dots). Values 
referring to a significancy level of 1, 2 and 3 a levels are reported 
too (dashed, dot-dashed and dot lines respectively). 



4. Next, the previous three steps have been repeated for a 
sample of 10^ mock lists of lines built in the following way. 
In order to take into account the varying number density of 
detectable lines along the Lyman-a forests, due to the vary- 
ing SNR, each forest has been simulated in chunks of about 
200 A. In each mock chunk, the number of simulated lines 
has been determined from a Poissonian distribution centred 
on the number of observed lines in that chunk, while the 
positions of the mock lines have been randomly generated 
following a uniform distribution within the related wave- 
length range of each chunk. The redshift intervals masked 
in the observed spectra were masked also in the simulated 
ones. The EWs of the mock lines have been randomly cho- 
sen among all the EWs measured by the fit of the lines in 
the observed spectra. In this way it has been possible, for 
each velocity bin, to compute the mean and the standard 
deviation of the number of occurrences for synthetic lists of 
lines. 

5. Finally, we have defined the three point probability excess 
(PES) as a function of the velocity difference, A v, according 
to the following formula: 



PE(v) = iV,,,(A^;)/iV,,^(A^;) - 1 



(7) 



The resulting PES is reported in Fig. [7| together with 
the 1, 2 and S a confidence levels. The PES is non-zero at a 
2 a level up to a velocity difference of ^ 250 km s~^. Most 
of the signal of the PES is due to the large number (26) of 
coincidences produced by the SS-S5-S6 QSOs triplet which 
is also the closest triplet (mean angular separation of 2.02 
arcmin corresponding to ^ 2 comoving Mpc). 

Fig. |8] shows the probability excess considering quadru- 
plets of Lyman-a lines. A significant signal at more than 
S a level is measured up to a velocity difference of ^ 250 
km s~^. Besides, one group of five coincident lines within 
100 km s-^ in the S2-SS-S4-S5-S6 QSOs is observed at a 
mean redshift of 1.825, an occurrence that has a probability 
P=0.01S to arise from a random distribution of lines. The 



portion of spectra where these five lines fall are reported in 

Fig.jg] 

In particular, it is possible to observe in the spectrum 
of the SS QSO the presence of a Damped Lyman-a system 
(DLA): the fitted Voigt profile gives a column density value 
of log A/" = 20.6. This DLA is associated with several metal- 
lic ion absorption lines found in the redder part of the spec- 
trum. Indeed at the same redshift we have found evidence 
of Civ, Fell, Siiv, Sim, Sin, Aim and Alii. This corre- 
lated H I absorption across lines of sight separated by ^ 14 

comoving Mpc could be interpreted as due to a gas fil- 
ament, but other explanations are possible such as detected 
Hi halos of clustered galaxies. 



4.2 The C IV lines coincidences 

We have also looked for correlation across two lines of sight 
for C IV absorbers. Given the characteristics of our sample, 
which allow us to investigate comoving scales larger than 
^ lh~^ comoving Mpc, we do not expect Civ coincidences 
due to the same C IV cloud, whose typical size should be of 
the order of tens of kpc (Ranch et al. 2001). Our lines of 



sight could pierce C IV gas connected with different galaxies 
in groups or clusters ( Francis et al.||2001 ). 

The procedure to identify CiV coincidences is similar 
but simpler than that previously employed for the Lyman-a 
lines. The first two steps are the same, with the difference 
that here lines closer than 500 km s~^ have been merged 
together. Then, we searched for coincidences of C IV merged 
absorbers lying along the redshif t space in windows smaller 



than 500 km s ^ . According to IScannapieco et al. ( 2006 ) 



this velocity scale corresponds to the characteristic size of 
coherent C IV structures. 

In order to establish the significance of our observa- 
tion, we have applied our algorithm for coincidences to 10^ 
groups of synthetic list of lines. The method for creating the 
synthetic list of lines is the same as adopted before in the 
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Figure 10. Portions of spectra representing the C IV systems found to be in coincidence between the S3 and S6 QSOs (left plots) and 
the Tl and T2 QSOs (right plots). The single CiV components are marked with increasing number. The same number is related to both 
the C IV lines. 



case of the mock list of Lyman-a lines. With the total set of 
simulated results, we could find the number N(n ^ no) of 
cases in which, in a random distribution of C IV absorbers 
along the lines of sight, the number of chance coincidences 
n is equal to or greater than the number of coincidences 
observed in our sample, no- Then the related probability P 
was simply computed using the ratio between N{n ^ no) 
and the total number of simulations, that is 10^. Finally, we 
computed the CL of our observation simply with the formula 
CL=l-P(n^no). 

We have found only two coincidences, in the T1-T2 and 
S3-S6 pairs, at a redshift of 2.062 and 1.486 respectively 
and separated by a comoving distance of 4.4 and 2.8 h~^ 
comoving Mpc. The statistical significance of these detec- 
tions are 90% and 85% respectively. Coincidences have also 
been found between the T2 QSO redshift and a T3 absorber 
(with a separation of 9.6 h~^ comoving Mpc) and between 
the S6 QSO redshift and a S3 absorber (with a separation 
of 2.8 h~^ comoving Mpc). These coincidences are in agree- 
ment with the fact that QSOs reside in over-dense regions of 
the Universe (e.g., [P'Odorico et al.||2008| ). No coincidences 
among absorbers have been found considering three or more 
lines of sight. The portions of spectra showing the two coin- 
cidences are reported in Fig. [lO] 



5 COMMON UNDER-DENSE REGIONS 

The two groups of QSOs in our sample are well suited to 
search for large under-dense regions in the IGM. Previous 
attempts to detect under-dense regions in the Lyman-a for- 
est have searched for regions without absorption lines along 
a single line of sight. They looked for portions of spectra 
where the observed number of absorption lines was signif- 
icantly smaller then the expected number from synthetic 
spectra ([Carswell Rees|1987||Crotts|1987]|Ostriker, Bajt- 



ence of an under-dense region along several close lines of 
sight. This would allow the detection of smaller under-dense 
regions maintaining a high confidence level. As done in the 
previous section, we focused our attention on searching for 
common under-dense regions among the lines of sight of the 
Triplet and among all the possible groups of lines of sight 
provided by the Sextet. 



lik Duncan|1988||Cristiani et al.|1995| |Kim et al.||2001 ) 

A strong improvement of the statistical significance of 
each detection can be provided by the simultaneous pres- 



Following |Rollinde et al.| ( |2003| , we can define an under- 
dense region along a single line of sight as a portion of the 
spectrum between and Ae , where the normalized flux / is 
higher than the specific threshold (/ — a/), where / is the 
average flux in the Lyman-a forest and erf is the value of the 
flux uncertainty. Note also that ai z ^ 2 regions above the 
mean density correspond on average to regions below the 
mean flux level of high resolution spectra ( Viel et al.|2008|). 
Since the SNR of the spectra analyzed by [Rollinde et ah] 
([2003j) is much larger than that of our spectra (15 — 70 vs. 
4—15 respectively), this simple procedure for the detection 
of under-dense regions cannot be adopted in the present 
case. Indeed, the fluctuations due to the flux noise cause 
the splitting of single large under-dense regions into several 
smaller regions, preventing their detection. 

In order to reveal the presence of under-dense regions 
in our spectra, we have developed the following method. 
The spectra are smoothed with a median filter character- 
ized by a very short kernel length, only three pixels. Af- 
ter the smoothing, common under-dense regions are looked 
for with the method described above (always adopting the 
same average flux and flux uncertainty). Then, the spectra 
are smoothed again and again searched for common voids. 
This process is iterated several times. The results are shown 
in Fig. [12] where we have reported the size of the bigger 
common under-dense region found in the Triplet and in the 
S3-S5-S6 group of Sextet QSOs as a function of the number 
of times the spectra have been smoothed. We report only 
the result concerning the S3-S5-S6 group because it is the 
only one in which a significant under-dense region has been 
detected. 
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Figure 11. Spectra of the two close groups of quasars (S3-S5-S6 on the left and the Triplet on the right) centred on the common 
under-dense regions. The extension of the common under-dense regions are shown with dashed vertical lines. 
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Figure 9. Portions of the spectra of the QSOs S2, S3, S4, S5 and 
S6 centred at 3435 A showing the aligned Lyman-o; lines. 
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Figure 12. Size of the common under-dense regions of Triplet 
group (empty triangle) and the S3-S5-S6 group (filled square) 
as a function of the number n of times the spectra have been 
smoothed. 



We can see how, after an initial linear growth of the size 
of the common under-dense region, a step behaviour occurs 
in both the series of data. Steps occur whenever two separate 
under-dense regions are merged together, that is when the 
absorption line (both real or spurious) lying between them 
is smoothed out by the n-th action of the median filter. We 
can notice also how, after ^10 smoothings, the size of the 
under-dense regions in both the Triplet and the S3-S5-S6 
group reaches an "asymptotic" value. Each asymptotic rate 
corresponds to the size of a common under-dense region en- 
closed between two strong lines (not necessarily belonging to 
the same spectrum) requiring a large number of smoothings 
to be cancelled. 
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A series of simulations has been carried out in order to 
evaluate the minimum value of the column density A/'(Hi) 
of the Lyman-a lines surviving after 6 and 8 smoothing pro- 
cesses in the Triplet and the S3-S5-S6 group, respectively. 
First, a mock spectrum has been created placing several lines 
characterized by an increasing logA/'(Hi), from 12 to 15 by 
steps of 0.05, and central wavelengths separated by 10 A. 
Line profiles have been modeled with a Gaussian, since we 
were interested only in the centre of the lines and not in 
reproducing correctly their wings. The Doppler par ameter h 
was set to a value of 20 km s~^(Kim et al. 2001); the de- 



pendence of the simulation results from this parameter is 
negligible. 

White noise has been added to the spectrum in order 
to obtain SNR — 7, roughly consistent with the observed 
SNR values of the Triplet, and of the S3, S5 and S6 spectra 
in the segments where the common under-dense regions are 
situated. Then, the length of the void has been measured, 
namely the portion of mock spectrum where / ^ (/ — a/), af- 
ter each smoothing. Since the lines are ordered by increasing 
column density and separated by 10 A, each increase in void 
length of 10 A corresponds to the vanishing of a specific 
line with known value of A/'(Hi). Repeating this procedure 
10^ times, we have computed the values, and the related un- 
certainties, of the column density of the lines washed away 
after a given number of smoothings. The effect of having dif- 



ferent SNR in the QSOs spectra (Fig. 11 ) instead of a fixed 
value of SNR = 7 is smaller than the error bars reported in 
Fig. [13] 

The same kind of simulations have been used to com- 
pare the action of the short scale filter iterated n times with 
respect to that of a median filter of length (2n+l) pixel em- 
ployed once. Both the series of data are shown in Fig. [13] 
with the related error bar, as a function of n. The trend 
followed by the two series at large n is different: for the 
short scale filter iterated n times the smoothing effect is 
relatively lower and the growth of the column density of 
the lines smoothed away is flatter. Furthermore, the related 
uncertainties are much smaller. This behaviour helps us to 
estimate more accurately the lower column density of the 
surviving lines, because the range of log NiYl i) between 12.5 
and 13 is sampled more accurately with respect to the upper 
curve. 



Fig. 13 (filled circles) shows that the asymptotic size of 
the voids observed in the Triplet and in the S3-S5-S6 group 
are referred to spectra where the absorption lines character- 
ized by logA^(Hi) < 12.8 have been smoothed away. 

The two largest under-dense regions common to the S3, 
S5 and S6 lines of sight and to the Triplet lines of sight are 
12.8 A (corresponding to 10.7 comoving Mpc) and 10.6 
A (or 8.8 comoving Mpc) in size, respectively. They 
are shown in Fig. [TT] The sizes of the related under-dense 
regions along the single lines of sight are: 13.4, 19.8 and 15.2 
A for the S3, S5 and S6 lines of sight; 10.6, 12.4 and 14.2 
A for the lines of sight of the Triplet. As already said, no 
common under-dense regions have been found considering 
any other combination of lines of sight in the Sextet. 

Assuming that the under-dense regions have a spheri- 
cal geometry (again following Rollinde et al.|2003 ), a void is 
then defined as the largest sphere that may be included in- 
side one connected region. We searched for the largest spher- 
ical region included in the three individual under-dense re- 
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Figure 13. Minimum values of the Hi column density of the 
Lyman-o; lines left after a given number n of smoothing with 
a three pixel long filter (circle) and after one smoothing with a 
(2n+l) pixel wide filter (square, slightly shifted along the x-axis 
for clarity). 



gions. This is equivalent to the requirement that the surface 
of the under-dense region matches the six edges of the three 
segments. We allow an uncertainty of few A in the position of 
the surface along each line of sight. This takes into account 
the fact that observations are in redshift space, so that the 
border of the under-dense region may have a peculiar ve- 
locity that shifts its position in wavelength, and that a real 
under-dense region probably does not have a simple bound- 
ary surface. This set of equations can constrain a sphere that 
has four degrees of freedom. Inside the sphere, the IGM will 
be under-dense along each line of sight and on its surface 
it will be over-dense at six positions identifying the edges 
of the individual under-dense regions. When we applied this 
procedure to the S3-S5-S6 group, the under-dense region can 
be parametrized with a sphere of 6.75 comoving Mpc 
radius. On the other hand, due to the positions of the mem- 
bers of Triplet in the sky, a spherical parametrisation of the 
void is not possible in this case. 

In order to assign a significance level to the detection of 
these two under-dense regions, we performed 10^ simulations 
of mock spectra reproducing the two triplets of QSO spectra. 
First, the lists of central wavelengths were created following 
the same procedure as defined in Section. 4.1. Then, the line 
profiles were created with the observed column density and 
Doppler parameter distribution (e.g., [Kim et al.,2002| ). Ran- 
dom white noise was added, according to the observed SNR 
of each spectrum in each chunk, and then the spectra were 
smoothed 8 times by the median filter. We verified that the 
results of the simulations do not change significantly for a 
number of iterations larger than six. Finally, we computed 
the distribution function of the probability of having one 
common under-dense region of given dimension along three 
lines of sight. Computing the random probability of having 
a common under-dense regions of size equal to or greater 
than ours, we obtained that the final significance levels of 
our detections are 91% and 86% for the Triplet and for the 
S3-S5-S6 group respectively. These two values are smaller 
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than the significance levels of the common under-dense re- 



gion detected by iRollinde et al. ( 2003 ) . We expect that this 



is due to the fact that we are dealing with only three lines 
of sight and that our spectra have a lower SNR. 



6 CONCLUSIONS 

In this paper, we exploited the capabilities of high resolu- 
tion UVES spectra of QSO pairs to study the 3-dimensional 
distribution properties of baryonic matter in the IGM as 
traced by the transmitted flux in the QSO Lyman-a forests. 
Our sample is formed by 21 QSO pairs evenly distributed 
between angular separations of 1 and 14 arcmin, with 
Lyman-a forests at a median redshift z ^ 1.8. By calcu- 
lating the correlation functions we compared the observed 
sample with a set of mock spectra drawn from a cosmo- 
logical hydro-simulation run in a box of 120 comoving 
Mpc, adopting the cosmological parameters of the concor- 
dance model. The simulated sample reproduces 50 different 
realizations of the observed sample (see Section 2 for de- 
tails). Furthermore, particular emphasis has been given to 
the search for alignments and other particular features of 
the Lyman-a forests in the two QSO groups present in our 
sample. 

In the following, we summarize our main results: 
Two-point statistics 

(i) The computed correlation functions are in substan- 
tial agreement with those obtained in our previous paper 
(D'Odorico et al. 2006" Paper I). There is consistency be- 



tween the clustering properties of matter in the IGM calcu- 
lated in the direction parallel and transverse to the line of 
sight using the parameters of the concordance cosmology to 
map the angular distance into velocity separation. This is 
also due to the relatively large error bars of the computed 
quantities. As an implication, peculiar velocities in the ab- 
sorbing gas are likely to be smaller than ^ 100 km s~^. 
Matter in the IGM is clustered on scales smaller than ^ 300 
km s~^ or about 4 comoving Mpc. The simulated corre- 
lation functions are consistent with the observed analogous 
quantities at the 1 — 2 a level for this particular sample. 

(ii) Thanks to the increased SNR for some of the spec- 
tra in our sample the enhanced clustering signal measured 
in Paper I with the cross correlation coefficient at a trans- 
verse velocity separation Av± ^ 500 km s~^ is no longer 
significant. 

Three or more point statistics 

(iii) Significant coincidences of Lyman-a absorptions 
have been detected among the lines of sight forming the 
Sextet implying the presence of coherent gas structures ex- 
tending ~ 14h~^ comoving Mpc. In particular an excess of 
triplets and quadruplets of lines within Av — 100 km s~^ has 
been measured at a significance of 16 and 9 a. Besides, one 
group of five coincident lines in the S2-S3-S4-S5-S6 QSOs is 
observed, an occurrence that has a probability P=0.013 to 
arise from a random distribution of lines. 

(iv) A method for the detection of under-dense regions 
in relatively low SNR spectra has been developed. One cos- 
mic common under-dense region has been detected in each 
QSO group; in the Sextet the under-dense region has a di- 
mension of 10.7 comoving Mpc and a significance level 



of 91%, while in the Triplet it has a dimension of 8.8 
comoving Mpc and a significance level of 86%. These values 
are significantly smaller than those typical of under-dense 
regions detected along single lines of sight. The under-dense 
region common to the lines of sight of the Sextet can be 
parametrized by a sphere of radius 6.75 comoving Mpc. 

This study of the cosmic web environment at z ~ 2 will 
soon be extended with many more QSO spectra either at 
medium resolution with the X-Shooter spectrograph (Kaper| 



et al.|2009 ) o r/and with the R ^ 10^ QSO spectra provided 



by SDSS-III dSchlegel et"aL]|2007D . The ultimate goal is the 
characterization of the topology of the IGM in real space in a 
variety of environments by using Lyman-a and metal lines. 
On this basis, both hydrodynamical simulations of struc- 
ture formation and high-resolution spectroscopic samples, 
like the one presented here, will provide the link between ob- 
servational quantities and the underlying density field and 
shed light on the impact of systematic effects. 
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